arXiv:nucl-th/9808035vl 13 Aug 1998 


Calculation of Higher Order Effects in Electron-Positron Pair 
Production in Relativistic Heavy Ion Collisions 


Kai Hencken and Dirk Trautmann 

Institut fiir theoretische Physik der Universitdt Basel, Klingelbergstrasse 82, 4056 Basel, 

Switzerland 


Gerhard Baur 

Institut fiir Kernphysik (Theorie), Forschungszentrum Jiilich, 52425 Jiilieh, Germany 

(August 13, 1998) 


Abstract 


We present a calculation of higher order effects for the impact parameter 
dependent probability for single and multiple electron-positron pairs in (pe¬ 
ripheral) relativistic heavy ion collisions. Also total cross sections are given for 
SPS and RHIC energies. We make use of the expression derived recently by 
several groups where the summation of all higher orders can be done analyti¬ 
cally in the high energy limit. An astonishing result is that the cross section, 
that is, integrating over all impact parameters, is found to be identical to the 
lowest-order Born result for symmetric collisions. For the probability itself 
on the other hand we find rather large effects at small impact parameters 
compared to the lowest order results, which translate to large effects for the 
cross section for multiple pair production. 
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Electron-Positron pair production in (peripheral) relativistic heavy ion collisions has 
attracted some interest recently due to the observation that at the relativistic heavy ion 
collider RHIC and LHC, the probability for this process calculated in perturbation theory 
violates unitarity, that is, gets larger than one, even for impact parameters of the order of the 
Compton wave length Xq ~ 386 fm. This fact was hrst shown in |^,^. The unitarity violation 
was then studied in a number of articles, taking into account higher order processes in pj-Q . 
It was found that the inclusion of these higher-order processes leads to the restoration 
of unitarity but also to new effects, mainly the production of multiple pairs. All studies 
also found that the probability for A^-pair production can be approximately described by a 
Poisson distribution. Therefore the probability from perturbation theory (in the following 
called “reduced probability”) has to be interpreted as the average number of pairs produced 
in one collision. Deviations from the Poisson distribution where studied for small impact 
parameters in P] and found to be rather small. 

Calculations of the impact parameter dependent probability in the external held approx¬ 
imation where calculated exactly in lowest order for small impact parameters b [0 and later 
also for all impact parameters The total cross section for pair production (that is, 

integrated over impact parameters) was also calculated. Using a Monte Carlo approach it 
was studied in [^. An analytical form of the differential cross section was found in [Q. 
Cross sections for multiple pair production were also given there. 

One open question is still the importance of higher order corrections — or even nonper- 
turbative effects — coming from the large value of the effective coupling constant Za ~ 0.6. 
Coupled channel calculations have been done at smaller energies (up to 7 ~ 2 in the collider 
frame). They always predicted a much larger probability compared to perturbation theory 
OQ. But recently the accuracy of these calculations has been questioned and by using 
a larger basis set smaller results were found. In the end results were only a factor of four 
larger than the perturbative results [T^, which is also in agreement with calculations using 


a spline-approach [Q . 

In a recent article the summation of the effect of the target to all orders was studied in 
the high energy limit (that is up to lowest order in I/ 7 ) for the related problem of bound-free 
pair production |^. Here the electron is created into a bound state of one of the ions. It 
was found that the summation to all orders could be done analytically and a fairly simple 
expression was found. The calculated probability for bound-free pair production was found 
to be slightly smaller than the hrst order calculations. 

A number of authors have in the mean time extended this approach also to the calculation 
of (free) electron-positron pairs. In a hrst article it was shown that here the summation 
can be done analytically again, leading to a rather simple modihcation of the matrix element. 
The authors of [O come to the same conclusion; they show also by integrating over the 


impact parameter that the cross section becomes identical to the lowest order Born result. 
In |]T3 the same conclusion is given. In the scattering of electrons in the held of colliding 


nuclei is studied, a problem which is closely related to pair-production. 

Of course there remain a few questions that still need to be addressed. The calculations 
were done using the Dirac-sea picture (that is, starting with an electron with negative energy 
in the initial state) and are therefore one-particle approximations of the full problem. The 
applicability of this approximation for situations where strong helds produce many particles 
has still to be investigated. The use of the reduced probability in the Poisson distribution to 
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get multiple pair production was only derived in the (usual) Feynman picture. Also only a 
part of all possible diagrams are included in the approach. Fig. |l|(A) shows typical diagrams 
which are included, whereas those of Fig. 0(B) are not. In |]^ it is argued that this second 
class of diagrams vanish in the limit 'j —>■ oo for electron scattering, therefore the neglect 
of them seems to be justihed. But such a rigorous proof is still needed for pair production. 
Finally it is well known from calculations with the approach of Bethe, Maximom and Davis 
that Coulomb corrections exist and are not small at the energies considered [P!|,P0|. This 
seems to be in contrast to the observation that the cross section should be identical to the 
lowest-order one for ion collisions. The most likely explanation of their absence is, that they 
are of subleading order in I /7 and were therefore dropped. 

It is our main aim to present a numerical calculation of the higher order effects in these 
expressions given in these articles in order to study their magnitude and their practical 
implications. Whereas the cross section for one-pair production is dominated by the large 
impact parameters, multiple-pair production is more sensitive to small impact parameters. 
Therefore we try to see whether a measurement of multiple pairs can be used to look for 
these higher order effects. 


We make use of the expression Eq. (54) of [|^. In addition we also keep all effects of 
hnite 7 in the expression. Written in our notation the (reduced) differential probability is 


P(p+,p_, 6 ) = / (fAqP{p+,p_, Aq) exp{iA^) 


( 1 ) 


where and p_ are the momenta of positron and electron and the Fourier-transform of the 
probability P is given for symmetric collisions {Z = Za = Zb) by: 

^(P+,P-, Ag) = ^ J d'^q 
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( 2 ) 


with p = Za, q' = q + Aq, Wi = (l,0,0,/3), W 2 = (1,0,0,—/3) and m the electron mass. 
This expression is almost identical to the one from perturbation theory, as given in [Q; the 
only difference is the additional exponents 1 for the photon propagator. This fact was 
already observed in [|^,0^. Integrating over b leads to a delta-function 6{Aq). Then the 
photon propagators are just complex conjugate to each other and only the absolute value 
enters. Therefore the cross section calculations of ||T0|,|6[| are exact in the high energy limit. 

Here we want to study the effect of the higher orders for small impact parameters. The 
total (single pair) cross section is completely dominated by the large impact parameters, 
especially in the high energy limit. Stronger deviations are to be expected mainly for small 
b, especially if b gets smaller than the Compton wavelength Xc ~ 386/m. 
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The new expression has a complex exponent, which makes the expressions oscillatory. 
Therefore a direct Monte Carlo integration is not possible. We rewrite it into a form with 
only standard Feynman integrals. We start by applying the usual Feynman-trick to group 
a product of two denominator into a single one, integrating over an auxiliary parameter. 
This trick is normally used for integer exponents. But it is easy to see by looking at a 
derivation of this in terms of 5-functions (see, e.g., f^), that the same expression also hold 
for complex exponents. We use it here in the following form: 


w 




w 


,1 — 2 ' 


'^dw 


B{1 + iri,l — ir]) Jo [wC + {1 — w)D]'^ 


( 3 ) 


Rewriting both photon propagators, we get two auxiliary integrations: wa and wb- In 
the new denominator the factor ±iT] just cancel and the remaining expression is of the form 
of (the square of) a propagator. The remaining two-dimensional integral over q can then be 
done in the same way as discussed in . We integrate over all hnal states of the electron 
and positron using VEGAS . 

Due to the oscillatory behavior of the numerator at the boundaries, a direct numerical 
evaluation of the auxiliary integrals is not useful. We therefore expand P for each Aq in 
terms of polynomials of the form [wa{^ — wa)]^[wb{^ — wb)]’'- Terms up to fc, / = 5 have been 
included in a £t, but convergence is already found with smaller exponents. The integrals 
over these polynomials can be expressed now in terms of the 5-functions. Our approach has 
the advantage that the coefficients of the polynomial expansion are independent of rj, apart 
from the trivial dependence. Therefore results for arbitrary rj can be calculated with no 
extra effort. Fourier-transforming this expression now with respect to Aq gives us 5(6). The 
total cross section, that is, integrated over 6 , is given directly by a = (27r)^5(Ag = 0 ). 

We have calculated P{Aq) and 5(6) for both SPS (Pb-Pb collisions at 7 = 10) and RHIC 
(Au-Au collisions at 7 = 100) conditions. As a check of the correctness of our calculations, 
we can compare them for 7^0 with the perturbative results in P], where a different 
approach was used. We get perfect agreement between those two, giving us confidence in 
our procedure. 

Figure ^ shows the Fourier transform for 7 = 100 and for different 7 . The effect of the 
higher orders is quite large, making 5 smaller up to about 30%. The shape of the curve 
itself is changed only slightly. For Aq —> 0 all curves coincide with each other. This has to 
be the case as the total cross section is identical to the lowest order one. 

Figure shows the impact parameter dependent cross section da/db = 27r65(6) for 
different values of rj. A deviation is only seen for small impact parameters, where it is quite 
large. 

Making now use of the Poisson distribution we can calculate probabilities for multiple 
pair production. For 7 = 100 we get the results shown in Fig. The higher order processes 
reduce the multiple pair production probabilities, but the probability for two-pair production 
is still large. Integrating over 6 we get the total cross sections as given in table |. For the 
single-pair cross section we use the approach of 0Q. It is clearly seen that the cross section 


^This could be improved by taking into account the Poisson distribution and also a lower bound 
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for multiple pair production is sensitive to the higher order effects in both cases. Therefore 
measuring them seems to be a practical way to look for these effects. 

In conclusion, we have calculated total probabilities for one and multiple pair production 
including higher order effects using the expression of [|T^. This approach effectively sums 


up all higher order diagrams in the high energy limit. We have calculated cross sections 
for multiple pair production and have found both probabilities and cross sections to be 
substantially smaller for the full calculation compared to the perturbative one. Therefore a 
measurement of this cross section should allow to really see these higher order QED effects 
in an experiment. 

Higher order effects can only be seen if one uses observables which are sensitive to small 
impact parameters. Multiple pair production is one such possibility. Depending on the 
experimental situation one could also think of other ways, e.g., the production of other 
particles together with an pair. 

In this article we have concentrated on “global properties” of the pair production, that is, 
total probabilities and cross sections. Our main aim was to demonstrate that calculations 
are possible. In order to see whether experiments will be able to see these effects, more 
differential studies are needed. We will present these and also details of the calculations in 
an upcoming publication. 


for the impact parameters. As these effects are rather small, we have neglected them here. 
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FIGURES 



FIG. 1. Diagrams of the type (A) are included in the matrix element. Diagrams of type (B) 
are assumed to be sub dominant for large 7 . 



FIG. 2. AqP(Aq), the Fourier-transform of the total probability is shown for 7 = 100. The 
trivial = (Za)^ dependence is divided out. Shown are the results for different rj for the full 
calculation. 
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FIG. 3. Shown is the differential cross section da/db for 7 = 100 and for different values of r]. 
The trivial = {Za)^ dependence is divided out. At small impact parameters the cross section is 
reduced quite substantially. 
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FIG. 4. The impact parameter dependent probability to produce N pairs is shown for up to 
three pairs. Shown are results for the lowest order Born result (stars) and also the full calculation 
(circles) both for 7 = 100 and Au-Au collisions. 










TABLES 


N 

Born (barn) 

full (barn) 

7 = 10, Pb-Pb {Z = 82,7] = 0.59) 

1 

4.21k 

4.21k 

1 

123 

84.4 

2 

8.61 

3.88 

3 

0.713 

0.212 

7 = 100, Au-Au (Z = 79,7] = 0.57) 

1 

34k 

34k 

2 

893 

624 

3 

113 

53.9 

4 

18.9 

6.04 


TABLE I. Total cross section for one and multiple pair production are given. Shown are the 
results for SPS/CERN and RHIC conditions. For the results for A = 1 the approach of ^ was 
used. 
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